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Coherent radiation spectrum from high energy e in a bent crystal with arbitrary curvature distri- 
bution along the longitudinal coordinate is evaluated, based on the stationary phase approximation. 
For a uniformly bent crystal a closed-form expression for the spectrum is derived. The spectrum 
features include a dip at its beginning and the sharp end, which may split into two breaks depending 
on the particle incidence angle. Estimates of non-dipole radiation and multiple scattering effects are 
given. The value for the crystal bending angle at which the dipole coherent bremsstrahlung theory 
holds best appears to be ~ 10 _4 rad. 

PACS numbers: 34.80.Pa, 29.27.-a, 41.60.-m, 78.70.-g 

Keywords: bent crystal, coherent gamma-radiation, multiple scattering 



I. INTRODUCTION 

Gamma-radiation emitted by electrons and positrons 
at their non-channeled passage through planarly oriented 
bent crystals has been investigated in a few recent exper- 
iments [l|, searching for signatures of the charged parti- 
cle volume reflection effect [3| in the radiation spectrum. 
However, the spectra observed were largely monotonous 
at typical photon frequencies u, and the difference be- 
tween the measured spectra from positrons and electrons 
was basically of order of the experimental errors. At first 
sight, that may appear surprising, since the inherent dy- 
namics and the final reflection angle of volume reflection 
itself are known to be sufficiently different in cases of 
positively and negatively charged particles. 

In article Q dedicated to computer simulation of the 
conditions of experiment [l[ it was mentioned that the ra- 
diation spectrum must contain a component of so-called 
coherent bremsstrahlung in a bent crystal (CBBC) (see 
Q, appendix), arising at fast charged particle highly 
over-barrier motion, when perturbative treatment of par- 
ticle interaction with the crystal is valid. Actually, radia- 
tion of that type may prove even dominant when the ac- 
tive crystallographic plane bending angle is many times 
larger than the critical value 9 C [2fJ - then the particle 
must spend most of its time traveling at angles to atomic 
planes much larger than critical, i. e., flying high above 
the potential barrier. Thereat, the frequency of the ra- 
diation emitted by the particle at a given instant is pro- 
portional to the local frequency of atomic plane crossing 
by the particle, as in ordinary coherent bremsstrahlung 
0. In course of the particle passage, the angle of atomic 
plane crossing varies, and the coherent radiation intensity 
accordingly sweeps over the spectrum. Close to the vol- 
ume reflection point, of course, the particle motion will 
become non-perturbative, but the radiation from that 
region contributes relatively little to the spectrum as a 
whole, in contrast to the situation with the elastic scat- 
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tering, where it is only the volume reflection point vicin- 
ity that contributes to the particle final deflection angle. 
The origin of the difference between the cases is obvious: 
the magnitude of the intra-crystalline transverse force is 
about the same all over the crystal, so all traversed crys- 
tal regions contribute commensurably to the total irra- 
diation energy. As for the particle deflection angle, it is 
sensitive to the force sign, and thus receives little contri- 
bution from the regions where the force oscillates rapidly, 
as it does at large angles of atomic plane crossing. 

In view of the described situation, prior to studies of 
radiation features stemming from non-perturbative seg- 
ments of particle motion in bent crystals, it seems expe- 
dient to determine the shape of perturbative CBBC spec- 
trum, which yields a wide continuous background and at 
the same time provides the conceptually simplest approx- 
imation. Unfortunately, so far it has not been evaluated 
in a form suitable for comparison with experiment. It is 
the purpose of this article to present a full, though basic 
calculation, and also to determine characteristic scales for 
the physics of the process. Last not least, we will assess 
robustness of the simplest CBBC theory against various 
deteriorating effects present in nature, such as multiple 
scattering and the dipole regime failure. It turns out that 
the range of the dipole CBBC theory is rather limited, 
although non- vanishing. 

In view of the universality of the CBBC radiation in 
bent crystals of various shapes (existing examples include 
cylindrically bent crystals, sine-shaped bent crystals 0], 
other microfabricated configurations may appear in fu- 
ture), we extend our treatment to the case of arbitrary 
crystal bend profile. Our ability to cope with it grounds 
on the applicability of stationary phase approximation 
allowing one to treat the crystal curvature as locally con- 
stant. That variant of the stationary phase approxima- 
tion is of different origin than the one arising in problems 
of synchrotron-like radiation, and does not contradict to 
the use of the dipole approximation in radiation. 

The paper is organized as follows. In Sec. |TT] we de- 
fine the bent crystal planar continuous potential and the 
corresponding transverse force. In Sec. IHII we proceed to 
evaluating the particle deflection angle by such a force 
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FIG. 1: Schematic of an ultra-high-energy particle passage 
through a thin bent crystal. The vicinity of point to of the 
trajectory tangency to crystal bent planes gives the main con- 
tribution to the particle deflection angle. 



to the leading order of high-energy perturbation theory 
of classical mechanics. In Sec. IIVI we evaluate the radia- 
tion spectrum in the dipole approximation, including the 
quantum effect of radiation recoil (allowing for photon 
energies to be of order of the initial electron's). In Sec.lVl 
the conditions of applicability of such an approximation 
are analyzed. Sec. I VII provides a summary 



z « —L/2 and z ~ L/2), and the particle essentially 
interacts with the continuous potential of the planes de- 
pending only on single coordinate x (see Fig. [T]). The 
distance between the bent planes is practically unaffected 
by the crystal curvature, and the equation defining each 
plane takes the form 

x p i(z) = Cpi + £(«), (1) 

constants C p i being equal-spaced with the inter-planar 
distance d. 

Then, if the continuous inter-planar potential in 
the bent crystal was Straight {%) (a periodic function 
with period d), which corresponds to an acting force 
Straight (x) = - aVet ** isbt , after bending of this crystal 
the force will modify to 

F(x,z)=F stia , igU (x-£(z)) (2) 

(still, it can be regarded as directed along x). For crystals 
of constant curvature 1221 . 



II. PLANAR CONTINUOUS POTENTIAL IN A 
WEAKLY BENT CRYSTAL 

A. Geometry definition 

At practice, for the coherent bremsstrahlung at over- 
barrier particle passage not to be spoiled by the particle 
multiple scattering on the target nuclei it is desirable to 
work with a crystal not thicker than a few millimeters 
(the same situation as for coherent bremsstrahlung in 
straight crystals). There are various techniques for man- 
ufacturing such crystals with bent atomic planes along 
the short direction, but basically they fall into two cat- 
egories. First - one of the crystal transverse dimensions 
is made short, and even shorter than longitudinal (i. e., 
<C 1 mm), which permits to bend the crystal along the 
longitudinal direction 0. Second - both transverse di- 
mensions are made sufficiently sizeable and the crystal 
so obtained is bent by some of the large dimensions, but 
securing that (short) crystallographic planes in the crys- 
tal, along which the beam is to be sent, acquire some 
bending, too. Specifically, the bending of the latter 
may be achieved through the action of anti-clastic forces, 
when the crystal is deformed simultaneously by both large 
transverse dimensions with different strength 
or one might just arrange the active planes to be under 
some angle a to the large face, and then they must bend 
along with the large face, although sin a times weaker. 
At the same time, for issues of particle passage through 
the crystal the deviations of the large crystal faces from 
planes may be neglected. 

In any case, for what concerns description of the par- 
ticle passage, the geometry implies particle incidence at 
some (small) angle 9q to z-axis, chosen normal to the 
crystal large faces (let the latter be located at positions 
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with R — const being the atomic plane bending radius. 
In what follows, we will rely on the stationary phase ap- 
proximation, in which the crystal curvature is treated 
locally, and is described by the local bending radius 



R(z) 
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to emerge naturally in the following. 



B. Nearly parabolic continuous potential and the 
corresponding force 



The dynamics of a high-energy particle in a crystal may 
be described by ultra- relativistic classical mechanics [ljj . 
As we have agreed, we will use perturbative description of 
particle interaction with the crystal; this is a rather com- 
mon approach in the theory of coherent bremsstrahlung. 
Conditions thereof will be specified later fSec. PVT) . 

In the perturbative treatment of classical particle pas- 
sage dynamics, as well as in quantum theory, it is advan- 
tageous to express the periodic continuous potential in a 
form of Fourier series. Such a representation is economic 
(provided only a few lowest harmonics dominate), and at 
the same time convenient when proceeding from descrip- 
tion of a straight crystal to a bent one. For evaluation 
of the particle trajectory and the emitted radiation, of 
direct relevance is not the potential but the force acting 
on the particle. To define the force - firstly, in a straight 
crystal - it is convenient to choose the origin of x-axis in 
the middle of some inter-plane interval, with respect to 
which the potential is an even function of x, whereby the 
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force has to be odd. Then, Fourier decomposition of the 
force involves sine functions only: 
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straight 



(x) 



^F n (-l)"sin 
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(the numerical factors have been introduced for further 
convenience). 

In the simplest but important case of (110) planar ori- 
entation of a crystal with diamond- type lattice (e. g., 
silicon), the inter-planar continuous potential is approx- 
imate by a parabola, and the corresponding force - by 
a linear-sawtooth function, whose Fourier decomposition 
reads 



F 



(110) 



(x) = — x 



n=l 



\x\<d/2 
(-1)" 



period. (4a 
2imx 



(4b) 



According to Pa)) . F± has the meaning of the force ex- 
tremal value achieved at x — > — | + 0, the sign of F\ being 
equal to that of the particle charge. On the contrary, in 
another important case of orientation (111), there are two 
different (but also nearly parabolic) wells within the pe- 
riod of the continuous potential [HJ , and the force Fourier 
decomposition turns somewhat more complicated: 
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Vt 



-V s 



(5) 



(Vl and Vs have meaning of depths of the alternate un- 
equal potential wells, while the well widths are exactly 
(Ll = jd and ds = \d). Anyway, once one factors out 
here the value of the first Fourier coefficient, 



T + Vs cos- 
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Eq. ([5]) will assume the form similar to (|4b 
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= 1 
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(7) 



Here c„ is a sequence of coefficients of order unity, neither 
increasing nor decreasing as n — > oo, and by definition 
ci = 1. [21 

To take into account thermal smearing of the poten- 
tial, i. e., the force continuity at the locations of atomic 
planes, the simplest though heuristic trick is to increase 
the power of n in the overall — factor of the trigonometric 
series: 




FIG. 2: Inter-planar continuous potential and the correspond- 
ing force (Eq. |(5J) shapes (in arbitrary units) for positively 
charged particles in a silicon crystals with (110) orientation. 
Solid lines: e = (cooled crystal, Eq. (|4a|) '): dashed lines: 
e — 0.4 (room-temperature crystal, Eq. l[8|l). For negatively 
charged particles the signs of the functions reverse. 



At that, the sequence c n (or its parameters Vl, Vs) may 
need to be corrected, but still the series is dominated by 
the first term, for which q ~ 1. Such a modification acts 
similarly to the conventional Debye- Waller exponential 
factor (which, in principle, is also heuristic, only its first 
order Maclaurin term being rigorously related to ther- 
mal averages) . We refrain here from discussing the exact 
relation of e with temperature T, only indicate that for 
the case of Si (110) at room temperature agreement with 
the potentials used in the literature is achieved at e w 0.4 
(see Fig. [2~]), whereas for Si (111) it takes e ~ 1. To supply 
more motivation to our ansatz, note that in what follows 
the summation of series of the type (JSJ with constant 
e and simple c n will yield Riemann (^-related functions. 
Such functions emerge as well for a zero-temperature po- 
tential (e = 0), only in the latter case the function argu- 
ments being integer or half-integer. Our approach cor- 
responds to extension of those arguments to arbitrary 
fractional values, i. e. to an "analytic continuation" , in 
order to model the effect of the temperature in a simplest 
way. None of the following numerical results (serving as 
estimates) depends crucially on this technique. 

Practical bent crystals are usually manufactured from 
silicon. The relevant physical parameters for silicon are 
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(9a) 
(9b) 



Note that product |-Fi|<i for those cases has practically 
identical values, which is important for the subsequent 
numerical estimates. But all our figures will be drawn 
only for a simpler case (110). 



III. INFINITESIMAL DEFLECTION ANGLE 



F, 



then 



= — L 

n=l 



(-l) n c n 2imx 
sin — - — , 



i.l+e 



e(T) > 0. 



(8) In this section we shall analyze the elastic scattering 

of ultra-high energy particles by the continuous force de- 
fined in the previous section. Choose the time reference 
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point at the moment of the particle passage through the 
middle of the crystal (z = 0), so that we may equate t ~ z 
(we will use units c = h = 1). Since the beam width in 
practice is always greater than the inter-planar distance, 
there is essentially a uniform distribution of particles in 
transverse impact parameters. Defining impact param- 
eter b of an individual particle as the trajectory initial 
asymptote intercept on x-axis, i. e. at z = (see Fig.[T]), 
the force acting on the particle in a bent crystal can be 
written as 




FIG. 3: Geometric interpretation of the coherence length of 
a fast particle traversing a bent crystal: Zefc is a half-chord 
within a curved crystalline plane (radius R) tangential to the 
next curved plane at distance d from the initial one. 



F(t) = e 



E 



2-"' 
(-1)" 



2F\ 
sin 27m 



b + 6 t-t(t) 



(10) 



with Q(v) - the Heavyside unit step function (zero for 
negative arguments and unity for positive ones). 

First of all, let us evaluate the particle deflection an- 
gle. Asymptotically, to leading order in the potential to 
energy ratio V/E, the deflection angle is proportional to 
the integral of force ([ID]) along the particle unperturbed 
straight path [24j |: 



quadratic terms, one brings (1111) to the form 
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Now, if point to belongs to the interval — j < t < -|, 
the integral converges in a small vicinity of this point of 
the width 
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where 



L/2 
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7ri? c ^ n 1+e 
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sin I 27rn- 

-L/2 



%t - m 



I* — *0| ~ ^EFC 

(t ) = ^2R(t (9 Q ))d, (15) 



where 



R(t) 



(16) 



dt, (11) 

(for a geometric construction see Fig. [3j) . Physically, this 
is the length on which the deflecting external field acts 
periodically, and will be referred to as the external field 
coherence (EFC) length. Then, the integration limits in 
(fT4f may as well be extended to infinity, and the integral 
(12) evaluates by a standard formula 



is the Tsyganov critical radius [13( (the above definition 
of R c yet allows it, along with F\, to have different sign 
depending on the particle charge sign). If the crystal 
bending is macroscopic, in the sense that displacement £ 
of the planes is (generally) 3> d, the integrand is rapidly 
oscillatory. For evaluation of such an integral, one may 
employ the stationary phase approximation 14]. This 
requires, in the first place, finding stationary phase points 
tn at which 



ts(0 O ) 



dt 



(13) 



i. e., the points of tangency of a ray with slope #o to the 
family of bent crystalline planes. For a convex function £ 
such a point is unique - and for simplicity we will assume 
this to be the case, dubbing it t (see Fig. [T]). Then, 
expanding function £(t) in Taylor series about to up to 



iAt' 



dt 



gif sgnA 



/ — , (ZmA = 0) 



with sgn function defined as sgm4 = — 1(+1) if A < 
0(A > 0). If, on the contrary, to falls beyond the in- 
tegration interval, the integrand is everywhere rapidly 
oscillatory, and the result is small (yet there are contribu- 
tions from the end-points, inferior to those from station- 
ary phase points, but they will be neglected throughout 
for simplicity). With this accuracy, 



rn (0 O ,6) we 
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+27T7J 



b + 9 Q to(0o) - Hto(Oo)) 



(17) 



5 



Hence, the magnitude of deflection angles is determined 
by the crystalline plane curvature radius in the point of 
the trajectory tangency to the bent planes. 

As for functions t (8o), £ (io(#o)) appearing in (fTT|). 
they may be obtained explicitly only if is a suffi- 
ciently simple analytic function. For instance, in case of 
a crystal of constant curvature ([3|) they express through 
the only available parameter - the plane bending radius: 



= -2*o. * 
(R = const) 



R 
2 l 



So, at practice reversible analytic parameterizations of 
the bending profile are favored. 

For the case of orientation (110), when c„ = 1, the sum 
in the right-hand side of (|17[) can be expressed in terms 
of Hurwitz (generalized Riemann) ^-functions: 



■ RJyJ 2 Rd 




FIG. 4: Deflection angle of a positively charged particle in a 
bent Si (110) crystal, in units of y/ ^ td , vs. the impact param- 
eter variable /3 = 4 + const (see Eq. (|f 9|H . Solid line: e = 
(cooled crystal), dashed line: e = 0.4 (room-temperature crys- 
tal). 
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where 
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b 

7r ( 



= 



b + 6 o t o (0o)-Z(t o (9o)) 
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(18) 



(19) 



characterizes the impact parameter of an oblique trajec- 
tory relative to the bent planes in point z = to, and 
C(a,w) is the Hurwitz zeta- function (l5j . while braces 
{. . .}f in its second argument indicate the entry fractional 
part (ranging from to 1). For orientation (111) the re- 
sult can be expressed through Hurwitz zeta-functions in 
a similar manner, but it is somewhat more bulky and 
shall not be quoted herein. 

Function (fl~8"l) (shown in Fig. 2J is not particularly sen- 
sitive to the value of e, except around the fracture points. 
The latter ones are located at 



2' 2' 
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In those points, function (fTS)) (and therewith (fTT)) ) is ex- 
tremal and achieves the value 
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lorn (#0, b) 

sgn£" 



= -e 
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- 1, 



where C( a ) 
function. 



3 W2fl(to(g ))d 
2 7 ttR c 



a is the ordinary Riemann zeta- 



The noticeable asymmetry of function(s) in Fig. [4] can 
be traced to the phase shift j, arising within the sta- 
tionary phase approximation. The average value of the 
deflection angle over the impact parameters is strictly 
zero, as an average of a sum of sine functions over their 
full period. The physical reason behind that is the uni- 
form distribution of an unperturbed particle flow over the 
crystal, entailing equal influence of positive and negative 
forces on the entire beam. So, in the adopted first order 
of perturbation theory there is no signature of volume 
reflection of the beam. Non-zero, however, is the angu- 
lar spread acquired by the beam, whose measure is the 
deflection angle mean square: 



(fllorn) = -j / db9 2 (0 O ,b) 



2R(t o (0 Q ))d 



T — 



3+2c ' 



(21) 



The latter quantity will also play part in the treatment 
of radiation. 



IV. DIPOLE COHERENT BREMSSTRAHLUNG 



Having established in the previous section the descrip- 
tion of the particle passage through the crystal, we are in 
a position to calculate the accompanying radiation. As is 
typical for bremsstrahlung from relativistic particles, the 
radiation is concentrated within a cone of angles ~ 7 _1 
about the forward direction, and those small angles may 
be treated inclusively. Let oj stand for the photon fre- 
quency (energy). Within the dipole approximation, the 
spectrum of radiation integrated over emission angles, 
and averaged over the impact parameters b of particles 
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in the beam, expresses through the acting force as [6j, [lj 
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where 



E' =E-w, 



<7minM 



2EE' 



(23) 



(allowing for uj ~ i£, to account for quantum radiation 
recoil effects), and 



dte lqt F{t,9 ,b). 



(24) 



Thus, q is a frequency of the force acting on the particle, 
and (|24p is the Fourier transformation of the intra-crystal 
force. 

With F(t) given by Eq. ||10p. its Fourier transform is 
conveniently evaluated by decomposing the sine function 
into a pair of exponentials: 

F(Q,6o,b) = — ^ nl+2e 

n=l 

L/2 d te^sin(2 7m b + 9 ° t -W 
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(25) 



When averaging the square of ([25]) over the impact pa- 
rameters, we employ the identity 
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by virtue of which 
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(27) 



(note the absence of interference between the exponents 
after the averaging). 

Evaluation of each of the two oscillatory integrals in 
(|27p is carried out by the stationary phase approxima- 
tion, as in the previous section. The external field co- 
herence length is the same as (fT5|) . only the location of 
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FIG. 5: Local frequencies of plane crossing (thick oblique 
lines), and the frequencies of higher harmonics (thinner 
oblique lines), as functions of the particle longitudinal coordi- 
nate (time). Drawn for a specific case of crystal with constant 
curvature, when q(t) dependencies are linear. Vice versa, 
the construction may be used for determination of station- 
ary phase points for a given frequency (dashed lines). Thick 
intercepts on the g-axis schematically show a double-step dis- 
tribution of the force frequencies described by function |F(g)| 2 
(see also Fig. [9J:). 



the stationary-phase point t n ±, about which function £(t) 
has to be Taylor-expanded, now depends on 9q and q/n. 
The equations for stationary phase points read 



t n ±(q,e ): £'(t n ±)-6 =± 



qd 
2ixn 



(28) 



(In accord with (|27p. radiation from different stationary 
phase points does not interfere). Physically, £' — do rep- 
resents the angle between the beam and the crystalline 
planes in the stationary phase point. Thereby, Eq. (|28p 
may be viewed as a local coherent bremsstrahlung con- 
dition, in which the local frequency of the driving exter- 
nal force is proportional to the local frequency of crys- 
talline plane crossing by the particle, which in turn is 
proportional to the local angle of the trajectory inclina- 
tion to the planes (cf. @)- Ultimately, the approximate 
^-integration gives 



V2 dte<>™^^) 



L/2 



»(~\tn ± M\) R{tn± ^ 60))d - (29) 



Substitution of ([27 1 H ty to (|2"2"j) leads to the result for 
the radiation spectrum of coherent bremsstrahlung in a 
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bent crystal: 
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2££' 



©( g -|W(?,«0)| )fl(*n+(«,0o)) 



+© ( 3 - |tn-(g,^o)| ) -R(*n-Mo))j>- (30) 

By virtue of the power factor n -s\'ii , this sum is 
strongly dominated by the term n = 1, so the coherent 
bremsstrahlung spectrum shapes for crystal orientations 
(110) and (111) appear to be only marginally different, 
and the temperature effect on the coherent radiation is 
also small. 

In what follows we will concentrate on application of 
formula (|30[) to crystals of constant curvature. In that 

case the solution to Eq. (j2"5j) is t n ± = (j) Q ± J R (see 

Fig. [5]), and in (j30|) one may draw constant R out of the 
integral, which allows one to accomplish the integration 
in terms of elementary functions. Introducing parameters 



(31) 



(signifying the frequencies of active crystalline plane 
crossing at the entrance and at the exit from the crystal) 
and a function 
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FIG. 6: Coherent radiation spectra (neglecting multiple scat- 
tering and incoherent bremstrahlung) , for a fixed value of 



<C E and several values of |#o 



Solid line: 9o 
dotted: \0 O \ -- 



= 

i 

R ■ 



(cj_ and lj+ coincide); dashed: |#o| = 3^; 
The impact of temperature effects on radiation is negligible 
the type of crystal orientation ((110) or (111)) mainly affects 
the frequency and intensity scales. 
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FIG. 7: Same as Fig. [6] for 
energy, the radiation spectrum shape is already strongly in- 
fluenced by the quadratic factor E' 2 /E 2 in Eq. (|33[) . obscuring 
the spectrum features. 



(0 < v < 1) 

the expression for the radiation spectrum assumes the 
form 



d^CBBC _ e 2 F 2 RdE r 



duj 



7r 3 m 2 E 2 



00 2 ( 

n—1 ^ 



/ <7n 



\ nq~ ' E 



+9(nq_ + q min )Q(nq+ 



)D 



+9(-nq_-q min ) 



D 



gmin W_ 

nq + ' E 



nq + ' E 

<Zmm W 



n\q-\ E 



(33) 



The behavior of the spectrum for different incidence 
angles is illustrated in Figs. [SJ [7J Despite being com- 
posed of discontinuous 0-functions, in total (|33| is ev- 
erywhere continuous, as is conditioned by its initial in- 
tegral representation (|3TJ)) . Still, there are discontinuities 
in the derivative of (|33p manifesting themselves as sharp 
curve breaks ("ankle" -type). Beyond the first two, major 
breaks [25] (corresponding to n = 1) the spectrum effec- 
tively ends, and only contributions from higher harmon- 
ics remain. Those main "ankles" are located at photon 
energies [26[ 



1 



U3± 



J- -I- 1 

E ^ 2 7 ^| g± 



(34) 
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Note that when |0 O | is only slightly below then 
cj_ <C lo+ , and at u> < cu- the spectrum develops a sharp 
spike, although superimposed on the background of equal 
height. That condition corresponds to a trajectory nearly 
tangential to the crystalline planes at the entrance or exit 
from the crystal, and although this feature may be of ex- 
perimental utility, one should beware that our present 
stationary phase approximation neglecting end-point ef- 
fects, as well as the dipole approximation itself, are in 
substantial error there. 

There are other important features of the CBBC spec- 
trum concerning dependencies on the geometric parame- 
ters R, L and 0q: 

(i) It is natural that the differential cross-section of co- 

herent radiation is proportional to the square of 
the field strength and to the square of the coher- 
ence length (fT5)) . We accentuate that the coherence 
length in our problem is of external origin and in- 
dependent of u - that is a length on which the 
particle-crystal interaction may be regarded as pe- 
riodic. In contrast, the photon formation length 
[27j Z form = > q^ 1 depends on w. It sets 

the scale of resulting photon energies, correlating 
with spacings q^ 1 between the neighboring crys- 
talline planes measured along the particle trajec- 
tory at the entrance to and at the exit from the 
crystal. 

(ii) One can check that as 10 — > the limit of (130)) is 



dEcBBC 

du 



2e^ 
to 3tt 



7 2 (0i ora ) 



(35) 



with (0| orn ) given by Eq. (|2Tj) . Apparently, this 
value does not depend on the crystal thickness L. 

(iii) The total radiation energy emitted per one electron 



E, 



CBBC 



da) 



dE c 



BBC 



duj 



(36) 



expresses rather simply and quite differently in two 
limiting cases: when the photon recoil effects arc 
negligible, and when they are crucial. If the "mod- 
erately high energy" condition 27 2 g + <C E is met, 
then g m ; n an d the second argument of all the 

D-functions in (|33p may be put to zero. For this 
case, one finds 



-Ecbbc 



8e 2 2 F 2 ~ 



3tt 2 



j2 n 2+2e ■ 

n=l 



md R , R 
^ << 2^L~ 1 °L 



(37) 
(38) 



In contrast to the differential intensity, the to- 
tal emitted energy here is proportional not to the 
square of the coherence length but to the crystal 
thickness. Remarkably, it does not depend on R, 



nor O , an d just equals to the total energy of coher- 
ent bremsstrahlung radiation in a straight crystal 
of thickness L. 

If the opposite condition 2j 2 q + » E is realized, 
then the first argument of all the D-functions in 
(|33l) may be put to zero, giving 



-Ecbbc 



dE, 



CBBC 



duj 



>o Jo 



du 



E 



19e 2 F?Rd 

367T 3 ' 



\ " 



j2 A^ n 3+2e 



E' 2 
E 2 







1- 



3a; 2 \ 
AEE'J 



2R~ % 



md R o R 



(39) 



(40) 



So, in this limit even the total radiation energy does 
not depend on the crystal thickness. 



(iv) At 



9o|» 



2R 



(41) 



(a large incidence angle or the straight crystal 
limit), Eqs. (|3TJ [HJ yield g_ ~ -q + < 0, w_ ~ 
u> + . Then in (|33|) the term with (nq_ — (faiin) 
vanishes. The next term containing 0(ng_ + 
9min)0("-<Z+ — qmin) is non-zero only in relatively 
small intervals n\q~\ < g m in < nq + , yet the corre- 
sponding D-function has its first argument close to 
unity and thereby is small (cf. Eq. (|32al0 : 



D 



nq + E 



1 



nq+ 



1 



2EE' 



< 1. (42) 



Lastly, the term containing 0(— nq- — q m i n ) con- 
tributes on the entire interval g m i n < i. e., 



basically at lj < ui- 



but there is a 



5 + 4^ T ^|0 | 

valuable cancelation between the corresponding D- 
functions: 



nq + E 



8 nf w 



-1-2- 



n\q-Y E / 

<7min / J_ 

q+ 



n\q- 



n 2 \q- 



Q- 



2ttL 



2EE' J nqi Rd 



(43) 



2tt|0 o | 



Therewith, the radiation spectrum reduces to 



dE T e 2 F 2 d 2 E 12 



In 



ymil 



1- 



2q n 



00 

\,j \l „>+*• V" »\q-\ //->/- 1 2 ' 2EE' ) ' 



(44) 
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which complies with the coherent bremsstrahlung 
spectrum in a straight crystal [6( (note that the de- 
pendence on R drops out). However, due to the 
9q 2 dependence of ([4"4"]l . with the increase of \6q\ to 
reach (|4Tj) the radiation intensity attenuates. Be- 
sides that, at large incidence angles the continuous 
potential approximation may be invalidated. 



CONDITIONS OF APPLICABILITY 



Our framework in the preceding two sections had been 
developed by the principle of maximal theoretical sim- 
plicity. In Sec. nn in our infinitesimal description of the 
particle deflection in the crystal we appealed to the high 
value of the particle energy. In Sec. IIVI we yet adopted 
the dipole approximation to radiation emission, which, 
however, is known [ljj to break down at a sufficiently 
high energy. Therefore, we have to investigate whether 
these two approximations are mutually consistent under 
conditions of a real silicon crystal, and if yes, what is 
their compatibility domain. Yet, besides the continu- 
ous potential influence on the particle there exists inco- 
herent scattering on individual nuclei, which affects the 
particle deflection as well as radiation. After all, in a 
case R \R C \ the condition of infinitesimal deflection 
certainly fails in vicinity of the volume reflection point, 
and that may also affect the radiation spectrum in some 
frequency domain. The present, last section comprises 
estimates of all the mentioned effects. 



A. Validity of the infinitesimal deflection 
approximation 

The condition of validity of the straight passage ap- 
proximation is the smallness of the particle transverse 
displacement relative to the inter-planar distance. Based 
on Eq. (|10l) . let us evaluate the particle transverse dis- 
placement as a function of time: 



Ax(t) = f dt' f dt" F{t ' r) 

J-L/2 J-L/2 



E 



— E 

71—1 



(-l) n Cn 



1+6 



* ,, r* ,//./'„ b+e t"-at") 

dt' I dt" sin [ 2nn d 



L/2 



L/2 



(45) 



Changing here the order of integrations, and again ex- 
panding in Taylor series about point to, one con- 
verts the double integral in (I45|) to a single one which is 
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FIG. 8: The particle local transverse departure from the ini- 
tial straight trajectory Ax(t), expressed in units of /| F c/-R c 
(solid line), and the local deflection angle 0(t), in units of 
Iefc/Rc (dashed line), as functions of the distance to the 
point of tangency to crystalline bent planes (for some spe- 
cific impact parameter). Both curves are built for the room- 
temperature case e = 0.4, crystal orientation (110) and for 
P = -0.7 (see Eq. 1(15)1). 



of Fresnel type: 
Ax(t)R c 



I 2 

'efc 



E 

n=l 



(-l)"c„ 



t—tn 

'efc 

■L/2+tp 

'efc 



dr 



t-t 



'EFC 



3in{27m(/3- t 2 )} . 



(46) 



It behaves as shown in Fig. [5] (by solid line). At large 
negative 

Ax(t)R c 

'efc 

t-t 



I 2 

'EFC 



1 



L/2 + 1 



1 ^ (-l)"- 1 ^ / / (L/2 + tpf 



I 2 

'efc 



^\t-tof 



E 



i-iy 



sin< 2nn (3 



{t-tpf 
I 2 

'efc 



(47) 



A curious feature here is the weak linear drift at the ini- 
tial stage, visualized in Fig. [5] and represented in Eq. (|4"T|) 
by the term proportional to 1 + L * 2+t • ^ ma y m " 
terpreted as a beam refraction at the entrance to the 
bent crystal. The refraction angle sign depends on the 
impact parameter at the entrance. However, by the ab- 
solute magnitude this effect is small, and for our current 
estimates less relevant. 

The most relevant is the behavior of function (|4"6"|) at 
t > to where it grows linearly (which corresponds to a 
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motion along the scattering angle final asymptote), as 
Ax(t)R c 



EFC 



->1 



*o) — 2^ -./^ sm(27TO/3 



7r * — ' n 3 / 2+e 



1 °° 



(-1)' 



,2+n 



■ cos < 2nn (3 - 



(i/2 + *o) 



'EFC 



I 2 

'EFC 



.(-ir 



87T 3 (t-t ) 



2 



,3+e 



— sin<27rn[ /3 



{t-tpf 
I 2 

'efc 



(48) 



In fact, the acting force contribution builds up only be- 
fore the particle reaches the hnal asymptote. For an 
actual estimate of the transverse displacement up to 
that moment one may simply take the value of (|48j) at 
~ 1 . For reliability of the straight passage approx- 
imation, the corresponding transverse displacement Ax 
needs to be less than the inter-planar interval half-width: 

Ax{t + Zefc) < d/2. (49) 

With the use of Eq. condition (gSJ boils down to 

(for perturb, defl. angle) (50) 
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This is physically obvious, since at R > |i? c | non- 
perturbative effects such as planar channeling, or volume 
reflection in the bent crystal already become important. 

Yet, for applicability of the stationary phase approx- 
imation the necessary requirement is that the external 
field coherence length Iefc — V2Rd be small compared 
to the crystal half-thickness, i. e., 



V2M< L/2. 



(51) 



In what concerns applicability of the straight passage 
approximation to the radiation, the corresponding con- 
dition is milder than (14U1 150")) . CBBC stems from the os- 
cillatory part of the particle motion on the initial or the 
final asymptote, represented by last lines of Eqs. (|47|) 
and (j48|) . Denoting that oscillatory part of motion as 
var Ax(t), the condition for the straight passage approx- 
imation applicability to the description of radiation is 
var Aa;L_ t , ^ , 2 <C d/2. Substituting for var Ax(t) the 

last line of (|47|) or (l48l) . one arrives at the requirement 



(52) 



L 2 2^2 d 
R 2>> 



Rr. 



That condition is essentially a product of (|50p and (I5T1) . 
so when (I5ip holds very well, inequality (|50[) may be re- 
laxed. 



deflection angle compared to the typical radiation angle 
7 _1 . So, let us evaluate from Eq. (JTUJ) the local angle 
of deflection from the straight path. Using again the 
stationary phase approximation, one gets 



dt 



,F(t>) 



-L/2 



'EFC * 

R c 7T ^ 

n—1 



(- l)"c, 



t-tQ 

( EPC 



drum {2im (/3 - t 2 )) 



This is an ordinary Fresnel integral; it is observed to 
converge within the range ZefCj an d its asymptotic forms 
at \t — to I S> ^efc are [28|] 
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(53) 



and 



9(t) 



'Born 



:»1 



'EFC 



R c (t-t ) 2tt 2 ^ 



)2 + f 



cos(27m(/3 — 



(t-tpf 
I 2 

'EFC 



(54) 



Now, for validity in weakly bent crystals of the dipole 
approximation to radiation, we need smallness of the os- 
cillatory part var j9(t) at \t — to\ < L/2. Substituting in 
Eq. |J53]) or the last line of Eq. ([54]) \t - t \ L/2, and 
replacing the sum by its typical value l/v2, we get 



Var ^l| t - t0 |<V2 <<;1 



L 
R 



(55) 

Here 9y is a parameter similar to 9y = ^ of [13]. [29[ 
With the use of parameters © , numerically one finds 



0.65-10" 



(56) 



both for Si (110) and Si (111). It is worth emphasizing 
that in a bent crystal the validity of the dipole approxi- 
mation to radiation depends on 9y smallness in compar- 
ison not with the Lindhard critical angle 



B. Validity of the dipole approximation for 
radiation 



2\Rr 



(57) 



Secondly, we had employed dipole approximation for 
the radiation, which presumes smallness of the particle 



(dependent on the particle energy via R c ), but with the 
active crystallographic plane bending angle L/ R. 
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C. Influence of multiple scattering on particle 
deflection and on radiation 



1. Deflection 

The angle of particle deflection in the continuous po- 
tential field also competes with the (rms, plane) angle of 
multiple scattering. The latter has a square root depen- 
dence on the medium thickness traversed tl8f : 



n (58) 

For electrons and positrons in silicon [181 ] 

'mult ~ 0.13 mm. (e^ in Si) (59) 

For multiple scattering not to affect significantly the par- 
ticle deflection in the target, 9 m uit(L) must be less than 
the angle given by Eq. (|T7|) : 



L \FiWRd 
< - — , 



D. Incoherent bremsstrahlung background 

Still another issue is that the coherent radiation re- 
ceives a background from incoherent radiation acts. A 
standard way to estimate the incoherent bremsstrahlung 
intensity in a crystal is to take the radiation in an amor- 
phous target made of the same material: 



dE 



BH 



L 

T 



4 / oj\ oj 2 

z\ 1 -e) + T 



<~>{E-oj). (63) 



Here Lq is the radiation length, for silicon amounting [1£ 



Lq = 9.36 cm. 



(64) 



The w-dependence of (|63[) is mild, and as an estimate of 
dE-QYi/doj one may take its value at oj ~ 0. 

To compare with, the spectral intensity of the CBBC 
radiation at an average radiation frequency oj ~ 
(see Figs. El [7]) is about half of its maximal value (f35l) : 



dE, 



CBBC 



du 



which entails 



2e 2 F?Rd ^ 4 

3 7r 3 m 2 n 3+2e 
n— 1 



(65a) 
(65b) 



R \ m 



Si (110) 
Si (111) 



(60) 



(|0 O | < L/2R) 
Numerically, Eq. (|65b|) gives 



2. Radiation 

Concerning the coherent radiation at typical frequen- 
cies, again, condition (j6"0)l is not relevant. Instead, one 
is to compare 9 mu \t with angle var#(t) at \t — to| ~ L/2. 
Should we be interested in the radiation angular distribu- 
tion, # mu it had to be count on the whole crystal thickness 
L. However, if only the (angle-integral) radiation spec- 
trum is looked at, for absence of the multiple scattering 
influence on it, angle var# at \t — to] ~ L/2 should be 
large compared to the multiple scattering angle only on 
the length 'efc : 



var# 



\t-t \~L/2 



mult {Iefc) 



(61) 



l. e. 



Z EFC ,,n R 

- — < ti v —. 

'mult J-i 



For the active crystalline plane bending angle this implies 
1.7' 



1.5 



10-4/ i^V /4 J Si (110) 



R 



Si(nr 



(62) 



3?r 
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(66) 



As had been mentioned at the end of Sec. IIV1 the coher- 
ent bremsstrahlung spectral intensity is independent of 
the crystal thickness L. 

For CBBC radiation to manifest itself prominently, it 
must exceed the incoherent bremsstrahlung contribution: 



dEcBBC dEBH 



doj 



du 



(67) 



With (jnUESl), it appears that the ratio gLEcbbc /c^bh 
depends only on the ratio L/R, i. e. on the active plane 
bending angle, with the proportionality coefficient 



dE, 



CBBC 



dEsH 



(68) 



E. Radiation at volume reflection (small oj domain) 



We had mentioned in Sec. lVBl that CBBC mechanism 
may be responsible for the generation of a large part of 
the radiation spectrum even when condition (|50[) is vi- 
olated. In the latter case, the infinitesimal deflection 
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highly over- motion close to highly over- 
barrier motion (half-) channeling barrier motion 




—\Ft\ 1 

highly over- motion close to highly over- 
barrier motion (half-) channeling barrier motion 

|F(q)| 2 




FIG. 9: (a) - exemplary graph of time-dependence of the force 
acting in a ((110) oriented, zero-temperature) bent crystal on 
a (positively) charged particle around the volume reflection 
point, at R S> R c (the figure corresponds to R — 50i? c ). 
The force discontinuities correspond to the particle passage 
through (sharp) potential maxima at the atomic plane po- 
sitions. In vicinity of point t = t vc n the trajectory draws 
nearly tangential to the maximum potential ridge, and in that 
sense is close to (half-) channeling, (b) - the same for nega- 
tively charged particles. The effective potential maximum is 
regular (F — > on top) and is situated approximately mid- 
way the atomic planes, (c) - schematic of the force Fourier 
transform modulus square. The dominant contribution to 
|F(q)| 2 and therethrough to dE co ^/duJ comes from the inter- 
val (jr v .r. < g < <?+■ 



approximation fails for evaluation of the particle final 
deflection angle, overestimating it, and hence the CBBC 
formula (|33[) must overestimate the radiation spectrum 
at sufficiently small uj, where it is proportional to the 
final deflection angle squared. Let us now estimate the 
scale of u! at which modification of CBBC radiation is 
needed. 

At R 3> |i? c | the actual mechanism of particle deflec- 
tion is volume reflection Q , whereat the magnitude of the 
deflection angle is of the order of Lindhard's critical angle 
(1571) . The contributing q- frequencies of particle oscilla- 
tion during the volume reflection are effectively bounded 
from below [301 ] by the value equal to twice the channeling 



dE coh /dw 

-r 2 <0L„>K- 
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FIG. 10: Schematic of a turnover in the coherent radiation 
spectrum due to the volume reflection. 



frequency 1/r: 



with 



T = T+ = 



q > q v .i-. = -, 

T 



(69) 



T+- 



V \Rc\d 

hlT-^T 



2tt 



(pos. char, part.) (70a) 
(neg. char, part.) (70b) 



(see Fig. IHt). The mode void below frequency arises 
because in the vicinity of the volume reflection point 
t = < rc fl (the closest approach to the axis of the crys- 
tal bending) the particle moves in each interval nearly 
by the channeling half period of the maximal amplitude, 
crossing thejxrtential ridges at a nearly grazing angle (see 
Figs. [9k, b) 31]. As a consequence, at radiation frequency 



1 



E 



(71) 



2rV 



the spectral intensity dE co h/duj related with |F((7)| 2 
through Eq. (|22|) . must have a turnover (see Fig. H0| . 



and drop at u -> Oto frT^.r. (with 0* r < (8l oln ), 
see Eq. ([71]) below) [l6J. Hence, at w ~ u v . T , there forms 
a spectral maximum, or rather a "hump" feature, since 
at uj > cj v r . the CBBC spectrum decreases rather slowly. 

A broad maximum similar to the one shown in Fig. [10] 
had first been discovered in computer simulations [2j, |4| 
of radiation at volume reflection, but its interpretation 
was not quite transparent. Now, we may conclude that 
as relative to CBBC, the volume reflection effect on ra- 
diation is only of suppressive, not enhancing character. 
It stems from the particle inability to sustain in a strong 
inter-crystalline field a quasi-periodic motion at too low 
frequencies - the over-barrier particle can not spend in 
an inter-planar channel a time longer than the channeling 
period (actually, half period). 

For the perturbative CBBC theory to have a significant 
applicability domain, frequency w v . r . must be much lower 
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than the CBBC spectrum end-point w + , which implies 

7ri m 



q v .t. < 



Rd' 27' 



(72) 



Condition g v r . <C ^ with numerical values ([9]) translates 



into requirement 



lOTeV, 



(73) 



which is guaranteed by the present accelerator capabil- 
ities, whereas condition g v r . <C basically, coincides 
with ([5"2"j) . The latter may also be regarded as a condi- 
tion for the particle energy (see Eq. (152)) below), which 
is, however, not very demanding. So, under the condi- 
tions of dipole radiation volume reflection effects should 
manifest themselves in a minor region of the spectrum, 
indeed. 

It is yet relevant to estimate the relative depth of 
the volume reflection dip. The actual value of the vol- 
ume reflection angle (in crystal with orientation (110)) is 

=s 9 r for 



|#v.r. I ~ f^c for positively charged, and \8 V 
negatively charged particles [lj|. Therefore, 



d£v.r.(0) 



dE, 



CBBC 



(~ 0) 



(vr/2) 2 
1 
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neg. char, part 
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2.5 



\Rc\ 
R 



(74) 



Thus, for the dip to be well discernible, one actually needs 
R to be at least a few times larger than R c . 



F. Crystal and beam optimal parameters 

Let us now assemble conditions (fSTI [521 15"5l 1^21 [67]) and 

examine their mutual compatibility, the variable param- 
eters being R and L. Eqs. ([55| 162]) together read 



tt V 'EFC 



which imply 



V^mult/^EFC ^ 1- 



(75) 



(76) 



This is reminiscent of the Landau-Pomeranchuk-Migdal 
(LPM) condition 



Ion 



(77) 



form 



< 



but at typical radiation frequencies one estimates I 
2g^ 1 ~ 2^, and according to inequality (|5"Tj) . it holds 
that ^form "C ^efc ■ So, the LPM condition appears to be 
less crucial than (|76|) . 

For fulfilment of condition (|76[) . with Z mu it fixed, one 
needs to have Zefcj i- e - -R and d, as low as possible. 
Note that the value of d is lower for orientation (110) 
than for (111), thus orientation (110) is more beneficial. 
But as for R, at practice it is normally at least in the 



range of decimeters, which gives \/Zmuit/^EFC — 4, while 
at highest R ~ 10 m one has v^muit/^EFC — 1-4. Thus, 
unfortunately, it is impossible to demand inequality (|76j) 
as really strong. Anyway, the optimal value for the active 
crystallographic plane bending angle is about 



L 
R 



EFC 



1/4 



1.3- 10" 4 . 



(optimal) (78) 



Then, the parameter of radiation non-dipoleness (also 
known as p-parameter [T?l |) is 



7 6»| 



-L/2 



e v ^- ~ 0.5, 



(79) 



whereas the parameter of radiation decoherence due to 
multiple scattering is about the same: 



A0 mu it(Z 



EFC J 



R 



\t-t \~L/2 



6 V L V I 



EFC 



0.5. 



(80) 



In view of the narrowness of condition (|76[) . it seems 
reasonable to suggest that since the size of the coher- 
ence length remains the same for all locations within a 
uniformly bent crystal, and the length Z mu it is constant 
as well, then even if condition (1521 fails (multiple scat- 
tering effects are substantial), the spectrum shape may 
still be roughly described by the present theory, only the 
intensity being suppressed by a factor depending on the 
ratio ^EFcAmuit- However, evaluation of such a factor is 
beyond the scope of the present article. 

Within our framework, presuming all the abovcmcn- 
tioned conditions to be fulfilled, let us check the last cru- 
cial condition (|oT)l, With ([75|. ratio (JBSJ will amount 
dEcBBC /dEBH ~ 7, which is satisfactorily high. 

Other relevant conditions ([5~T1 [521 are rather easy to 
fulfil. At bending angle (|T8[) Eq. (|5ip demands for the 
crystal thickness 



- > — w 6-10 
d L 



12 Urn i Si(110) 
(81) 

Eq. (|52|) together with (f78|) sets the lower bound for the 
electron energy: 



lGeV. 



(82) 



Towards experimental investigation of CBBC itself, let 
us yet determine the full set of parameters at which 
CBBC effects are least deteriorated and the sharpest. 
As we had mentioned, lower R are favored for cleaner 
CBBC, but R is tightly interrelated with L, which must 
satisfy inequality (|81[) . Taking, marginally 



50 
80 



Si (110) 
Si (111) 



(83) 



(such thin samples are available 
entails through Eq. ((78j) 



R* 



cf. Rcf. 1 of 10]), 



0.4 
0.6 



Si (110) 
Si (111) 



(84) 
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For what concerns the electron energy, to avoid a pro- 
nounced volume reflection turnover one needs, according 
to ([74| . smaller ratio R/\R C \, i. e. higher E. Besides 
that, if one wants the "ankle" features in the radiation 
spectrum not to be smeared out by the photon recoil ef- 
fects (cf. Figs. |6]and [7]), one should arrange the condition 
< E, which under ([75} implies 

„ m 2 dR f 300 1 „ , r f Si (110)1 , oc , 
S «^L^{ 500 j GeV - { Si (111)} ( 85 ) 

If in marginal fulfilment of ([55} one takes 

£*«150GeV, (86) 

we derive 

whereby R+/\R C \+ « 1.6, and according to Eq. ([741 . the 
dip should not develop. Parameters ([531 [541 156")) are 
"ideal" for checking the calculations in the present pa- 
per; the spectrum thereat must look like that in Fig. ([S} . 

The oj range is up to ; , 1 nd ■ ~ 35-^50 GeV and the ex- 
's + 2^L 

pected spectral intensity maximum is dEc ^, c ^ s» 0.04. 
The most demanding condition seems to be the initial 
electron beam collimation degree S9q < jji ~ 70 //rad. If 
not achieved, an averaging of the radiation spectrum over 
the electron beam incidence angles must be performed. 

VI. SUMMARY AND DISCUSSION 

The present study substantiates the notion that spec- 
tral intensity of radiation from ultra-high-energy elec- 
trons and positrons in a bent crystal is a sum of con- 
tributions from particle motion intervals on which the 
local angle between the particle velocity and the bent 
crystalline planes is definite, and so each such contribu- 
tion is similar to coherent bremsstrahlung in a straight 
crystal. The length of an elementary coherence interval 
is ~ /efc — V2i?d, implying that (i) the crystal must be 
much thicker than the mentioned value (see Eq. ([51]) ); (ii) 
the radiation spectral intensity (Eqs. (|30H33I) ') is propor- 
tional to igpc °c Rd, and does not depend on the crystal 



thickness. Still, under the "moderately high energy" con- 
dition ([3"51) the crystal thickness determines the spectrum 
extent, and therethrough, the total energy emitted. 

The characteristic feature of CBBC from a single elec- 
tron is the well-defined end of the radiation spectrum, 
whose position depends on the active crystallographic 
plane bending angle L/R. At an angle of electron inci- 
dence on the crystal comparable to jfe this end of spec- 
trum splits into a pair of breaks (see Fig. [5}. That feature 
must in principle be experimentally verifiable with a suffi- 
ciently well collimated initial beam (59q < jji)- The best 
experimental conditions for investigating CBBC were de- 
scribed in Sec. IV Fl 

We have also qualitatively discussed the modification 
of the coherent radiation spectrum in the domain of small 
u> owing to the onset of the volume reflection phenomenon 
possible when R 3> |i? c |. That modification is of purely 
suppressive character and manifests itself as a dip at the 
beginning of the spectrum. Next to the dip, around fre- 
quency (|7ip there appears to be a maximum in the spec- 
trum, but it is not to be interpreted as a resonance. 

The theoretical description adopted in this article had 
resorted to many simplifications - it did not properly in- 
corporate the temperature dependence of the potential, 
neglected multiple scattering, and relied on an infinitesi- 
mal approximation to the particle deflection (in the bulk 
of the medium) as well as on dipole description of the 
radiation. Conditions ([ST1 [781182} under which those ap- 
proximations hold, altogether appear to be restrictive for 
the crystal bending angle (see Eq. ([78} ). so generalization 
to a non-dipole treatment, and an account of the multi- 
ple scattering would be highly desirable. Nonetheless, let 
us mention that the dipole CBBC conditions are quite 
nicely met, e. g., in recent experiment Q. Comparison 
of the CBBC theory with the available experimental data 
is intended elsewhere. 

In conclusion, let us remark that although our paper 
presumed dependence of the crystal deformation only 
on one, longitudinal, coordinate, in principle higher- 
dimensional deformation cases are conceivable, emerging 
under application of torsion, or owing to intrinsic crys- 
tal mosaicity. In those cases the stationary phase ap- 
proximation must still be applicable, but the description 
should inevitably become more sophisticated. 
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